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Abstract 

Starting from the kinetic equations for the fluctuations and correlations of a dilute gas of inelastic 
hard spheres or disks, a Boltzmann-Langevin equation for the one-particle distribution function 
of the homogeneous cooling state is constructed. This equation is the linear Boltzmann equation 
with a fluctuating white noise term. Balance equations for the fluctuating hydrodynamic fields are 
derived. New fluctuating forces appear as compared with the elastic limit. The particular case of 
the transverse velocity field is investigated in detail. Its fluctuations can be described by means 
of a Langevin equation, but exhibiting two main differences with the Landau-Lifshitz theory: the 
noise is not white, and its second moment is not determined by the shear viscosity. This shows 
that the fluctuation-dissipation relations for molecular fluids do not straightforwardly carry over 
to inelastic gases. The theoretical predictions are shown to be in good agreement with molecular 
dynamics simulation results. 

PACS numbers: 45.70.Mg, 45.70.Qj, 47.20.Ky 
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I. INTRODUCTION 



The (modified) nonlinear Boltzmann equation for the one particle distribution function 
provides an accurate description of transport phenomena in a low density gas of inelastic 
lard spheres or disks . These particles are often used to model granular fluids 

6], specially in the rapid flow regime [?J]. The Boltzmann equation does not provide any 
direct information about correlations and fluctuations in the gas, other that the particle 
velocity moments. Nevertheless, methods used in the derivation of the Boltzmann equation 
have been extended to obtain kinetic equations for the equal and different time correlations, 
in the same low density approximation. The general idea is that in order to obtain these 
equations the needed approximations are the same as those used to derive the Boltzmann 
equation itself. 

One of the earliest and physically more transparent methods to study fluctuations is 
;hat of Langevin equations. Almost 40 years ago, in a seminal paper Bixon and Zwanzig 
8j showed how a Boltzmann-Langevin equation could be constructed by generalizing the 
reasonings leading to the Boltzmann equation for molecular gases. The latter describes 
the behavior of the average value of the one-particle distribution function, while the former 
incorporates the effects of the fluctuations. As the authors indicated themselves in the paper, 
the derivation was based on physical intuition and analogy. A more systematic derivation 
of the same result, starting from first principles, was given in ref. |9|. 

A second approach to the study of correlations in dilute gases makes use of functional 
analysis. Its more general result is a kinetic equation for a generating functional at low 
density, from which all multi-point correlations can be obtained by functional differentiation 



11 



10 ] . A closely related general scheme for the study of correlations is the hierarchical method 



121 ]. The starting point are hierarchies of coupled equations for the time distribution 



functions describing the fluctuations and correlations. Then, the hierarchies are closed by 
using the same kind of approximations as needed to derive the kinetic equation, i.e. the 
Boltzmann equation in the case of dilute gases. This method has been recently extended 
to describe fluctuations and correlations of dilute inelastic gases in their simplest state, 
the homogeneous cooling state (HCS) ]13]. As an application, the fluctuations of the total 
energy were studied and a good agreement between theory and simulation results was found 
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One of the aims of this paper is to translate the above formalism in terms of kinetic 
equations for the correlation functions into a Langevin equation formulation, i.e. to extend 
the fluctuating Boltzmann equation to the case of inelastic hard spheres or disks. The 
relationship between kinetic equations and the fluctuating Boltzmann equation has been 
analyzed in detail in molecular gases [3, Q]- One advantage of the Langevin formulation 
is that it is closer to the fluctuating hydrodynamic equations. Actually, the fluctuating 



171 ] to lead to the same 



Boltzmann equation for molecular systems has been shown 

ay 

Langevin equations for the hydrodynamic fields as obtained by Landau and Lifshitz 
using thermodynamic fluctuation theory. The noise terms in these equations are assumed 
to be white with second moments determined by the Navier-Stokes transport coefficients 
of the fluid. Their expressions are known as fluctuation-dissipation relations of the second 
kind |l9]. 

The derivation of fluctuating hydrodynamic equations from the fluctuating Boltzmann 
equation for inelastic hard particles, will be also addressed here. Attention will be focussed on 
a particular state, the HCS, and on a specific hydrodynamic field, the transverse component 
of the velocity. The main conclusion will be that the fluctuation-dissipation relation for 
elastic gases can not be directly extrapolated to inelastic ones, but it needs to be significantly 
modified. The second moment of the noise is not determined by the Navier-Stokes shear 
viscosity. Moreover, the noise can not be assumed to be white. These theoretical predictions 
are in qualitative and quantitative agreement with molecular dynamics simulation results. 

The consideration of the HCS does not imply by itself that the results obtained here are 
not relevant for other states more accesible experimentally. The HCS plays for inelastic gases 
a role similar to the equilibrium state for molecular gases. In the case of molecular systems, 
the expressions of the transport coefficients obtained by linearizing around equilibrium are 
the same as those appearing in the nonlinear Navier-Stokes equations as predicted by the 
Chapman-Enskog method and successfully used in many far from equilibrium problems 20]. 
Also, the fluctuation-dissipation relations derived for near-local-equilibrium states in the 



origma 

states 



Landau and Lifshitz theory have proven to be accurate for many other hydrodynamic 



21] . For dilute gases composed of inelastic hard particles, the equivalence between 



the transport coefficients obtained by linear perturbations of the HCS and by applying the 
Chapman-Enskog procedure has also been established [22] . Something similar might be 
expected for the fluctuations and correlations. 



3 



In the system being considered here, the particles move freely and independently between 
consecutive collisions. More specifically, they are not coupled to any external energy source 
or thermal bath, contrary to the driven granular gas models. For these models, the linear 
resp onse to an external perturbation [23( as well as the validity of the Einstein relation 
24( have been investigated by numerical simulations, and some empirical models have been 
proposed. It is not evident a direct relation between the free model considered here and the 
above driven models. 

The plan of the paper is as follows. In Sec. [Til t ne kinetic equations for the one-time 
and two-time correlation functions of a dilute gas in the HCS derived in ref. jl3] are shortly 
reviewed. These equations are translated into an equivalent Boltzmann-Langevin equation 
for the one particle distribution function in Sec. IHIl When written in the appropriate 
variables, this equation is the linear Boltzmann equation to which a fluctuating force term 
is added, similarly to what happens in molecular elastic gases. An expression for the second 
moment of the fluctuating force in terms of the collisional Boltzmann kernel is derived. In 
Sec. IIVI the fluctuating hydrodynamic fields are defined, and balance equations for them are 
obtained from the Boltzmann-Langevin equation. They involve formal expressions for the 
fluctuating pressure tensor, the fluctuating heat flux, and the fluctuating cooling rate. In 
addition, an intrinsically inelastic fluctuating force shows up in the equation for the energy. 

To get a closed description for the hydrodynamic fluctuations, expressions for the heat 
flux, the pressure tensor, and the cooling rate in terms of the fluctuating hydrodynamic 
fields are needed. This can be accomplished by means of the Chapman-Enskog procedure. 
Here only the case of the transverse component of the velocity field will be considered. As 
a consequence, only the expression for the non-diagonal elements of the pressure tensor 
is required. This is computed in Sec. |VJ The final result is a Langevin equation, that 
is the linear macroscopic equation for the transverse velocity field plus a fluctuating force 
term. Therefore, the structure is similar to what one could expect by extrapolating from 
the corresponding equation for molecular systems 25| . Nevertheless, the noise term is not 
white and its second moment is not given by the usual fluctuation-dissipation relation. It 
is verified that the obtained theoretical predictions are in good agreement with molecular 
dynamics simulation results. Section IVIII contains some general comments and conclusions. 
Finally, the appendixes provide some details of the calculations needed to derive the results 
presented in the bulk of the paper. 
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II. KINETIC EQUATIONS FOR THE HOMOGENEOUS COOLING STATE 



The system considered is a dilute gas of TV smooth inelastic hard spheres (d = 3) or disks 
(d = 2) of mass m and diameter a. The position and velocity of the ith particle at time t 
will be denote by Ri(t) and Vi(t), respectively. The effect of a collision between particles i 
and j is to instantaneously modify their velocities according to the rule 



V = V 



a 



{a ■ Vij) 



(1) 



j -j ■ 2 

where = V — Vj is the relative velocity, a is the unit vector pointing from the center of 
particle j to the center of particle % at contact, and a is the coefficient of normal restitution. 
It is defined in the interval < a < 1 and it will considered here as constant, independent of 
the velocities of the particles involved in the collision. A more realistic modeling of granular 
gases would require to consider a velocity dependent restitution coefficient 5]. 

Given a trajectory of the system, one-point and two-point microscopic densities in phase 
space at time t are defined by 



N 



F 1 (x u t) = ^ t S[x 1 -X j {t)] 

3=1 



and 



N N 



(2) 



(3) 



F 2 (x 1 ,x 2 , t) = J2J2 5 ^ Xl ~ 6 t X2 - ' 

i j^i 

respectively. Here = {Ri(t), V(t)}, while the Xj = {r^, v^} are field variables refer ring 



to the one-particle phase space (/x space). The density Fi(xi,t) obeys the equation [ill. 1 1 31] 



d 



d 

dri 



F 1 {x 1 ,t)= / dx 2 T(x 1 ,X2)F 2 (x 1 ,X2,t), 



(4) 



with 



T(xi, Xj) = a 



d-l 



da Q(vij ■ a) 



\vij ■ tr\ [a 2 5(r ij - a)bj-{v it Vj) - <J(r - + <r)] , (5) 

where da is the solid angle element for a = a /a, r 12 = V\ — r 2 , G is the Heaviside step 
function, and v 2 ) is an operator replacing all the functions of V\ and v 2 to its right 

by the same functions of the precollisional values vl and v\ given by 

bZ X ' 



Vn 



- Vl = vi — [a ■ v 12 )cr, 

la. 

,_i l + a.^ 

b„ v 2 = v 2 + — — (cr • v 12 )a. 



2a 



(6) 



It is seen that Eq. for F\ involves the two particle density F 2 . Actually, it is the first 
equation of an infinity hierarchy [ljj. 

The averages of Fi(x\,t) and F 2 {x\,x 2 ,t) over the initial probability distribution of the 
system p(T, 0), T = . . . ,X^}, are the usual one-particle and two-particle distribution 
functions, 

fi(xi,t) = (Fi(xi,t)), f 2 (xx,x 2 ,t) = (F 2 (xi,x 2 ,t)), (7) 

where the notation 

(G} = J dTG(T)p(T,0) (8) 

has been employed. Two-time reduced distribution functions can also defined from the 
microscopic densities and the initial probability distribution. The simplest one is the two- 
particle two-time distribution function, 

= {F 1 {x 1 ,t)F 1 {x' l ,t')). (9) 
From the definitions in Eqs. (j7|) and Q it follows that 

fi^Xijtjx'^t) = 5(xi -x[)fi(xi,t) + f 2 (x 1 ,x' 1 ,t). (10) 
It is convenient to introduce one-time and two-time correlation functions by 

g 2 (x 1 ,x 2 ,t) = h(xi,x 2 ,t) - /i(xi,t)/i(x2,t), (11) 

and 

/ii il (ari,t;a/ 1 ,t / ) = f hl (x u t; x[,t') - fxixu^f^x'^t'), (12) 
respectively. Equation ffTUl) translates into 

h lil (x 1 ,t;x[,t) = 8{x x - x' 1 )fi(x 1 ,t) + g 2 (xi,x[,t). (13) 

In the low density limit, a closed set of kinetic equations for fu 92, and hx t \ can be 
derived |l3| by extending the methods developed for molecular gases jl2j ]. They can be used 
to analyze the average properties as well as correlations and fluctuations in arbitrary states 
of a dilute granular gas. Here attention will be restricted to a particular state of a freely 



evolving granular gas, the so-called homogeneous cooling state (HCS) [26(. Macroscopically, 



it is characterized by a uniform number of particles density n, a vanishing velocity field, 
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and a uniform time-dependent temperature T(t). It is further defined by the one-particle 
distribution function having the scaled form [l| 



f(v,t) = nvv d (t) x (c), 



where 



v {t) 



2T(t) 



m 



1/2 



(14) 



(15) 



is a thermal velocity and xi c ) is an isotropic function of the scaled velocity c = v/vo(t). 
The distribution x(c) and the granular temperature T(t) are specified by the pair of coupled 
equations 



dT 



Cod_ 
2 dc 



In the above expressions, 

Co 



'1 — crW 2 



(cx) = Jc[c\x}- 



dci / dc 2 cl 2 x(c 1 )x(c 2 ) 



2T (^±5) d 

is the dimensionless cooling rate in the time scale s defined by 



(16) 
(17) 

(18) 



Jo 



A 



(19) 



with A = (na x ) , and J c [c|x] is the inelastic Boltzmann collision term. Its explicit form 



is 



Jc[c\x] = / dciT (c,Ci)x(c)x(ci), 



(20) 



T (c, c x ) = j dSQ[(c- Cl ) • a\{c - Cl ) ■ a [a^b' 1 ^ c x ) - l] . (21) 

The variable s defined in Eq. ( JT9l) is proportional to the accumulated number of collisions 
per particle. For thermal velocities, i.e values of c of the order of unity, a good approximation 
;o the solution of Eqs. ( TL61) and (TTT1) is provided by the first Sonine approximation, in which 
IB 



with 



x(c) = e ^j- 2 [l + a 2 (a)S^(c 2 )} 



n d/2 

c 4 d + 2 



4 



c 2 | d(d + 2) 



(22) 



(23) 
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and 



a 2 {a) 



16(l-a)(l-2a 2 



9 + 24d + (8d - Al)a + 30a 2 - 30a 3 



In the same approximation 



r (d/2) d 



3a 2 (a) 
16 



(24) 



(25) 



A numerically exact solution of Eqs. (130]) and ffTTj) has been recently reported in 



28|. The 



wo-particle one-time correlation function of the HCS is assumed to have also a scaled form 



£ 2 (n 2 , «i, v 2 , t) = nX \ 2d (t)g(l 12 , ci, c 2 ), (26) 

where the scaled length scale I = r/X has been introduced. The dimensionless correlation g 
does not depend on s and obeys the equation 

d 



C12 ■ 



01 



A( Cl ) - A(c 2 ) 



12 



^(iu, ci, c 2 ) = 5(Z 12 )T (ci, c 2 )x(ci)x(c 2 ) 



(27) 



where A(Cj) is the linearized Boltzmann collision operator 291 ] . 



A(cj 



Co ^ 

rfc3To(ci,c 3 )(l + P i3 )x(c 3 ) - — — • c { . 



(28) 



The operator Pij interchanges the labels of particles % and j of the quantities to its right. 



For the two-particle two-time correlation function the scaling reads 



13| 



hx,x{xx, t; x[,t') = nX-\ d {t)vo d {t')h{lx - l[, c 1; s - s'; c[) 



(29) 



and the kinetic equation is 

9 d 1 ~ 

- + ci • — A(ci) /i(Zi - l[, cx, s - s; c[) = 0, 

valid for s > s' > 0. The initial condition for this equation is 



(30) 



h(lx-l'x,cx,0;c[) = hx,x(h - l'x, cx; c[) 



g(lx - l[, ci, c[) + 6(cx - c'JSih - l[)x(ci] 



(31) 



An equation for this distribution follows from Eqs. (|T71) and f l27l) . 

Cl "dT 1 +c ' 1 'W 1 ~ A(ci) " A(c ' l} l{h ~ Vl ' Ci; c ' l} = 5(Zl ~ (ci ' c ' l); (32) 
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with 

f ( Cl , c'O = - [A(ci) + A(c' 1 )] (J( Cl - ci)x( Cl ) + T (c 1; c'Jx^Ox^). (33) 

Equations fl27|) and fl30l) describe the correlations between fluctuations in the HCS. They 
become closed once the solution to Eqs. (fTBT) and (|T7|) is known. In the next section, an 
alternative and consistent description to that provided by these kinetic equations will be 
developed. 



III. FLUCTUATING BOLTZMANN EQUATION AROUND THE HCS 

Equation (j4]) is an exact consequence of the dynamical equations governing the motion 
of the particles. The aim of this section is to approximate it in such a way that give a closed 
description of the effective dynamics of a dilute granular gas in the HCS. To do so, the 
spatial separation between the centers of colliding particles will be neglected in the operator 
T(xi,x 2 ), and F 2 (xi,x 2 ,t) will be approximated by an effective (Boltzmann) two-particle 
phase space density at the mesoscopic level F 2 B (xi,x 2 ,t). Moreover, the dimensionless time 
scale s and length scale I introduced in the previous section will be used. Then, Eq. (jlj) 
becomes 

d~s + J ^' Cl + Cl '^)^ l(Zl ' Cl ' s)= / dc ^o(ci,c 2 )Fi(h, Cl ,h,c 2 ,s), (34) 
where dimensionless phase space densities have been defined by 

Fxih, ci, s) = rr^WFifa, t), (35) 

F 2 B (h, Cl , l 2 , c 2 , s) = n-\ M (t)F 2 B (x 1; x 2 , t). (36) 
Comparison of the ensemble average of Eq. (I3~4l) with Eq. f[T7|) gives the conditions 

(F 1 (/ 1 ,c 1 , S )) H = x(ci), (37) 

J dc 2 T (c 1 ,c 2 )(F*(l 1 ,c 1 ,l 1 ,c 2 ,s)) u = J c [c\ X }. (38) 

The subindex H in the angular brackets indicates that the ensemble average is taken over 
the probability distribution for the HCS. 

The deviation of the microscopic density from its average value is defined by 

SF^h, c u s) = F 1 (l 1 , c x , s) - x(ci). (39) 
9 



An evolution equation for this quantity follows by subtracting Eqs. (134")) and (|T7|) . 



9 ^ 9 , 5 l c / ■ , 1 



dc 2 T (ci,c 2 ) 



F<f (Zx, d, h, c 2 , s) - x(ci)x(c 2 ) 



(40) 



The structure of this equations suggests to introduce a cluster decomposition for F 2 S of the 
form 



F 2 (h, ci, Zi, c 2 , s) = x(ci)x(c 2 ) + x(ci)£Fi(li, c 2 , s) + x(c 2 )£Fi(*i> Ci, s) + $ 2 (Zi, c 1; c 2 , s). 

(41) 

This equation defines the microscopic correlation density $f(Zi, Ci, c 2 , s). Substitution of 
its ensemble average in Eq. (|38|) yields 



dc 2 T (ci,c 2 )($ 2 3 (Zi,Ci,c 2 ,s))h = 0. 



(42) 



Moreover, use of Eq. (14~TI) into Eq. fj4*0l allows to rewrite the equation in the equivalent form 

5F 1 (l 1 ,c 1 ,3)=S{h,c 1 ,8), (43) 



9 9 A / \ 



<9s 



where 



S(h,cx,s) = I rfc 2 T (c 1 ,c 2 )$ 2 3 (c 1 ,c 2 ,Z 1 ,s) 



(44) 



and the operator A(ci) was defined in Eq. ff28l) . Equation (j43l can be interpreted as a 
fluctuating Boltzmann-Langevin equation for the one-particle distribution function g), 
171 ]. with the "noise term" being S^Zx, c 1; s). Of course, this does add any new physical 
insight by itself in the understanding of the starting equation (134j) . The relevance and 
usefulness of this representation will depend on the properties of the noise term. A first one 
follows directly from Eq. fj42|) . that is equivalent to 



(S(l 1 ,c 1 ,s)) H = 0, 



(45) 



i.e. the noise has zero average. In the following, other properties of S will be derived by 
requiring consistency with the results derived in the previous section. Multiplication of Eq. 
( )43|) by SFifti, c[, s') with s' < s, followed by averaging gives 



9 9 */ N 

^ + Cl .--A(c ly 



(5F 1 (h, ci, s)5F 1 (l[, dy, s')>h = {S(h, ci, ^^(Zi, c£, s')>h- (46) 
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From the definition of 5Fi(li, c 1; s) it is easily verified that 

(5Fx(h, c h s)5F 1 (l[, c[, s'))h = n^X'^ih - l[, c h s - s'; <), (47) 

where h is defined in Eq. ( )29l) . Therefore, consistency of Eqs. (|46|) and ( 1301) implies that 

(S(l 1 ,c 1 ,s)5F 1 (l' 1 ,c l 1 ,s')) K = 0, (48) 

for s > s'. Since, by hypothesis, the parameters of the system are such that the HCS is 
stable, the long time solution of Eq. f|43l is 



5F 1 (l 1 ,c 1 ,s)= f dTe^ L( - h ^S(h, Cl ,T), (49) 

J — oo 

where the linear operator 

L^cO^AfcO-d- J- (50) 
has been introduced. Using Eq. ( ]49l . it is obtained 

[L(h, Cl ) + L(l[, c[)} (SFxih, Cl , s)5F 1 (l[, c[, s')) u 

= - f dTe^ L ^\S{l l) c l) r)S{l' 11 c! l1 s)) }1 

J — oo 

- f C /re (s - T)L(r -< ) (5(Z 1 ,c 1 , S )S(r i ,c' 1 ,r)) H . (51) 

J — oo 

This equation must be compared with Eq. fl32l) . having in mind Eq. f|47|) . The time inde- 
pendence of the right hand side of Eq. ( 1321) prompts to introduce the hypothesis that the 
noise term S is Markovian, and write 

(S(h, c 1; s)S(l[, c[, s'))h = #(ci, ci)5(i! - - s')- (52) 

On introduction of this into Eq. fl5T|) and comparison with Eq. fl32l) . it follows that 

H(c 1 ,c' 1 ) = n- 1 \- d T(c l ,c[), (53) 

with T defined in Eq. (1531) . 

The properties given by Eqs. ( l4"5l) . (HE]), and ( |52l) guarantee that the description provided 
by the Langevin equation ( 1431) leads to the same expressions for the two-particle, one-time 
and two-time correlation functions as the formulation in terms of reduced distributions 
functions reviewed in Sec. |TTJ 
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IV. FLUCTUATING HYDRODYNAMIC FIELDS AND BALANCE EQUATIONS 



The fluctuating number of particles density, N(r,t), momentum density, G(r,t), and 
energy density, E(r,i), are defined in terms of the microscopic phase space density as 

N(r,t) = j dv F 1 (x,t) } (54) 
G(r,t) = J dvmvFifat), (55) 

/771V ^ 
dv—F x {x,t). (56) 

Dimensionless deviations from their averages values in the HCS are given by 

5p(l,s) = = f dc5F l (l,c,s), (57) 



Sut(l, s) ee SG ^fl = I dcc5F 1 (l, c, s), (58) 
mnv {t) J 

Se (l,s) ee 25 = \ / dcc 2 SF 1 (l,c,s). (59) 



dnT{t) d 

The quantity 5u)(l, s) is the dimensionless velocity field. Balance equations for these fluctu- 
ating fields follow by taking velocity moments in the Langevin-Boltzmann equation ( f4*3i) and 
using the properties of the noise S. Some details of the calculations are given in appendix 
lAl The resulting equations read 

— 5p(l,s) + --5u(l,s) = 0, (60) 

(|-|) M i lfl) + |. ff n(, ia ) = o, (6i) 

^ - Co) Se(l, s) + ^ • s) + 5Co(i, s) + ~ • <ty(I, a) = S e (l, *)• (62) 

In the above equations, 8U(l, s) and 8(f)(l, s) are the fluctuating pressure tensor and heat flux, 
respectively. Their definitions in terms of the fluctuating one-particle distribution function 
are 

SU(l,s) = ^^\ + J dcA^SF^^s), (63) 
6<f>(l,s) = J ' dcX(c)6Fi(l,c,s), (64) 



12 



where I is the unit tensor of dimension d, and 



A(c) = cc— \, 



, 9 d 2 + 2 
£ c = c 2 — 



c. 



(65) 



(66) 



The term 5(o{l, s) represents the fluctuations of the cooling rate about its average value in 
the HCS. Its formal expressions is 



5( {l,s) = ^ ° ]^ / cfc! / dcac^xCci)^^,^, s). 



r(^)d 

Finally, S^Z, s) is a fluctuating force term having the properties 



(67) 



(S t (l,s)) H = 



(68) 



and 



{S e {l,s)S e {l',s')) n = 7C 1 \- d 6(s-if)6(l-l') 



—^2 Co I dcc^xic) 



dci / ^c 2 x(ci)x(c 2 )$(c 1 ,c 2 ) 



(69) 



with $(ci, c 2 ) given by Eq. (1A16[) . This noise term is intrinsic to the inelasticity of collisions 
and has no analogue in normal fluids. Of course, in the elastic limit a — *■ 1, xi c ) becomes 
a Gaussian and the fluctuating force S e is seen to vanish i n ag reement with the well known 



results for hydrodynamic fluctuations in molecular fluids [18j. The other main differences 
between Eq. (1621) and the one for molecular gases is the presence of the two terms involving 
the cooling rate, (o, and its fluctuations, 5(o- The presence of these terms is directly associ- 
ated with existence of the cooling term in the macroscopic equation for the average energy 



V. LANGEVIN EQUATION FOR THE VELOCITY FIELD 

Equations fl60|) - fl62|) do not provide a closed description of the hydrodynamic fluctuations 
of a granular gas in the HCS until 5U, 8(f), and 5(q are expressed in terms of the fluctuating 
hydrodynamic fields. This turns out to be not a simple task, and attention will be restricted 
in the following to the equation of the velocity field 5u)(l, s), Eq. flBTT) . 
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Given two functions f(c) and g(c), their scalar product is denned as 

(f\g} = J dc X -\c)f + (c)g{c), (70) 

where f + {c) is the complex conjugate of /(c). Next, a projection operator V is introduced 
by 

d+2 

Vf(c) = ^ p {e)^\f). (71) 

13=1 

Here, the functions £p(c),/3 = 1, .. . ,d + 2 are the eigenfunctions of the linear Boltzmann 
operator A(c) defined in Eq. ff28l) . corresponding to the hydrodynamic part of its spectrum. 
Therefore, they are solutions of the equation 

A(c)fr(c) = A^(c). (72) 



Their expressions are 22|, [29 ] 



ei(c)=x(c) + ^-Mc)], 6(<0 = -^, 6(c) = • [cx(c)\ ■ (73) 



The associated eigenvalues are found to be 



Co , Co 



Ai = 0, A 2 = |, A 3 = -|, (74) 
the eigenvalue A 2 being <i-fold degenerated. Finally the functions £g(c) are 

?i(c)=x(c), I 2 (c) = c X (c), Uc)=(^ + l)x(c). (75) 



The sets of functions {^(c)} and {£g( c )} are seen to have the biorthogonality property 

<^'> = <W'> (76) 

/3, j3' = 1, 2, . . . , d + 2. This guarantees that P as defined by Eq. (ITTil is really a projector 
operator, i.e. it verifies V 2 = V. It projects any function of c onto the subspace spanned by 
the hydrodynamic eigenfunctions of A. 

In the following, it will be more convenient to work in the Fourier representation. The 
Fourier transform of 8Fi(l, c, s) is 

<JFi(fe,c,s) = J dl e~ ik l 5Fi(Z,c,s). (77) 
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By means of V, SFi(k,c,s) can be decomposed into its hydrodynamic and non- 
hydrodynamic parts, 

SF^k, c, s) = V5F x {k, c, s) + V±5F 1 (k, c, s), (78) 

where T-^ = 1 — "P. The Fourier representation of the balance equation for the velocity 
fluctuations, Eq. ( 16T1) . is 

{~L ~ 2) 5UJ ^ S)+lk ' <5n(fc ' S) = °' (79) 

8U(k, s) = S > I + J dck{c)5F x {kiCis)i (80) 

where A(c) is defined in Eq. (1651) . Getting an explicit expression for Sn(k,s) in terms of 
the fluctuating hydrodynamic fields is the next issue to be addressed. By direct evaluation, 
it is easily verified that 

J dcA(c)&(c) = 0, (81) 
(3 = 1, . . . , d + 2. Hence Eq. (ISO)) is equivalent to 

5U(k, s) = MM> 1 + J dc A(c)P ± 5F 1 (fc, c, s). (82) 

To compute V±Fi(k, c, s), the operator Pi is applied to both sides of the Fourier transform 
of the Boltzmann-Langevin equation (|43|) . 

9 



P ± [A(c) - ik ■ c] } V±SF 1 (k, c, s) = -V±ik ■ cPSF^k, c, s) + V±S(k, c, s), (83) 
as J 

where use has been made of the property V±AV = 0, that is a consequence of the fact that 
P projects over a subspace generated by eigenf unctions of A. The solution of the above 
equation can be formally written as 

P±SF 1 (k, c, s) = U(k, c, s)V±SFx(k, c, 0) + f ds'U(k, c, s')V± \-ik ■ cVbF x {k, c,s- s') 

Jo L 

+S(k,c,s-s')\, (84) 

with 

U(k, c, s) = exp [sV±L(k, c)Pj , (85) 

L(k,c)=A(c)-ik-c. (86) 
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Taking again into account that the HCS is assumed to be stable for the system considered, 
the first term on the right hand side of Eq. can be neglected for large enough times 

s. Moreover, to derive hydrodynamic equations valid to Navier-Stokes order, only the first 
order in k of the pressure tensor is needed. To this order, 

/ ds'U(k, c, s')V± (-ik ■ c) VSF^k, c, s - s') 
Jo 

-> / ds'e s ' r ^ Aic)r ^V ± (-ik-c)V5F 1 (k 1 c 1 s-s') 
Jo 

= [ ds'V ± e s ' A ^(-ik-c)P5F 1 (k,c,s-s'). (87) 
Jo 

Then, for large s Eq. reduces to 

V±_8Fx(k,c,s) = [ ds'V ± e s ' A{c) (-ik- c)V8F x (k,c,s - s') 
Jo 

+ I ds'U(k,c,s')V±S(k,c,s- s'), 
Jo 



(88) 



and substitution of this into Eq. (!82|) yields 



where 



and 



5U(k, s) = 6 ^ k ' ^ I + ^n(fc, s) + R(fe, s), (89) 



5xn(fc,s) = J ds' j dc A(c)e s ' A ^(-tk ■ c)V5F x (k, c, s - s'), (90) 
R(k,s) = J ds' J dcA(c)U(k,c,s')V±S(k,c,s-s'). (91) 



Upon writing Eq. ( 1901) . Eq. ( IHTj) has been employed to remove the operator Vj_ appearing in 
the first term on the right hand side of Eq. (188|) . Because of the isotropy of the operator A(c), 
only the projection onto £2(0) gives a non- vanishing contribution to the above expression 
for 5 1 n(fc, s), that can be simplified to 

6 1 n(k,s) = f Q ds ' f ^A(c)e s ' A(c) (-^-c)^(c).^(c)|^ 1 (fe,c,s-s')) 
= J ds' j dcA(c)e s ' A{c \-ik ■ c)£ 2 (c) ■ 5u(k,s- s') 
~ J ds' J dcA(c)e s ' A{c) (-ik ■ c) e - s ' Co / 2 £ 2 (c) • 6u>(k,s). (92) 
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In the previous transformations, the definition in Eq. (!58|) has been used, and it has been 
taken into account that to lowest order in k, 



5cu(k,s-s') = e- s '^ 2 Su(k,s), 



(93) 



according to the balance equation for the fluctuating velocity field, Eq. flBTl) . Using again 
the symmetry of A(c), it is obtained: 



5ill ij (k, s 



-irj{s) 



ki5ujj{k, s) + kj5u>i(k, s) — — 5yk ■ 5u>(k, s) 



(94) 



This is the same as the Navier-Stokes expression for the pressure tensor with the only 
difference that the average macroscopic velocity field is substituted by the fluctuating one. 
It involves the (time-dependent) dimensionless shear viscosity rj(s) defined by 

1 



rj(s) 



d 2 + d-2 



E E f ds ' j dc ^(c)e^ A -^ 2>i (c) Cj . 



(95) 



This expression is equivalent to the one obtained from the nonlinear Boltzmann equation 
for inelastic hard spheres or disks by the Chapman-Enskog method 



30 



Green-Kubo formulas derived in ref. 



31| and also to the 



321 ]. Let us remark that the results obtained here 



apply in the limit of large s. It is in this limit when hydrodynamics in the usual sense is 
expected to apply. If this is true, the shear viscosity in Eq. (1951) will become independent 
of s. Although there is no a mathematical proof of this "ageing to hydrodynamics" for 
granular gases up to now, numerical evaluation of the right hand side of Eq. (I9~5l) by 
using the direct Monte Carlo simulation method has shown the existence of such a limit 
value 33J . Moreover, the simulation results for the shear viscosity rj are in good agreement 
with the expression obtained by evaluating the Chapman-Enskog result in the first Sonine 



approximation 



30 



3l|, 



77(a) = [8u(a) - Co(«)] 



-1 



7T 2 



via) 



-(3- 3a + 2d) (1 + a) 



1 - 



a 2 {a) 
32 



(96) 
(97) 



2y/2d(d + 2)r {d/2) 

When Eq. ( |94"|) is substituted into Eq. (|89|) and the resulting expression is used into Eq. 
9]), a Langevin-like equation is obtained for the velocity field, 



^-|) 5u>(k,s) + ~5e(k,s)k + ?j 



d — 2 

k 2 5u)(k, s) H — k ■ Suj(k, s)k 

Uj 



W{k,s), (98) 
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with the noise term given by 



W(k,s) = -ik ■ R(k,s), (99) 
where the term R(k, s) is defined in Eq. (l9Tj) . It follows from Eq. (j4"5]) that 

(W(M))h = 0. (100) 



A formal expression for the correlation function of W is obtained directly by using Eq. 
Its conversion into an explicit one, valid to Navier-Stokes order and, therefore, consistent 
with the left hand side of Eq. (|98|) . will be carried out in the next section for the particular 
case of the transverse component of the velocity field. 

VI. THE NOISE TERM IN THE EQUATION FOR THE TRANSVERSE VELOC- 
ITY FIELD 

The transverse part of the fluctuating velocity field, 5cv±(k, s), is defined by 

5oj ± (k, s) = <Wfe, s) - 5w(k, s)-—k. (101) 

k 

Its evolution equation can be written down directly from Eq. ( 1981 . 

^ - | + rfk^j 5u> ± (k, s) = W ± (k, s), (102) 
where 

W ± (k,s) = -ikh- — \ :R(k,s). (103) 
Substitution of the expression of R given in Eq. (l9~Ij) yields 



W±{k,s) = -i J ds' J dck-c(c-^^kjU{k,c,s')V±S{k,c,s-s'). (104) 

By using this expression, it is shown in appendix [B] that for two components, W± t i(k, s) and 
W±j(k, s), of the noise of the transverse velocity field, to lowest order in k it is 

(W ± , i (k,s)W ±jj (k',s')) H = 5 hJ 5 k ,- k ,^rk 2 G(\s - s'|), (105) 

for s > s' ^> 1. Here V is the volume of the system in the length scale I, i.e. V = N/n\ d , 
and 

G(\s\) = J dc J dc'A :E j / (c)A :cy (c')^HCs(c,s;c'), (106) 
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with 



*Ahcs(c, s; d) = J dlh(l,c,s;d). (107) 



The distribution h(l, c, s; d) is defined in Eq. (1291) . Then, by integration of Eq. (130]) it follows 
that ^hcs( c , s; d) obeys the equation 



9 \ / \ 

as - A(c) 



V) H cs(c,s;c') = 0, (108) 



valid for s > 0. In principle, this equation must be solved with the initial condition 
V^hcs(c; d), given by 

^hcs(c;c') = J dlg(l,c,d) + 5(c-d) X (c), (109) 



obtained by integration of Eq. (l31j) . Nevertheless, it has been shown in ref [33] by particle 
simulation that contributions from the correlations in the HCS of dilute granular gases are 
negligible, at least for not too strong inelasticity, a > 0.5. Therefore, the term involving g 
in Eq. (11091) has been neglected in the results reported below. 

The solution of Eq. ( 11021) in the limit of large s can be written as 



where 



_,i(M= f d Sl e x ^ s -^W ±!i (k,s), (110) 

J — oo 



A ± (fc) = |-^ 2 . (Ill) 



Next, using Eq. (11051) it is obtained that 

(5uj ±ii (k, s)5uj ±ti (k', s)) H = k 2 5 k ,- k ,yj-j, (112) 
with the coefficient rf defined as 

r? = 2 I" dsG{\s\)e£. (113) 
Jo 

This coefficient can be computed in the first Sonine approximation. Some details of the 
calculations are given in appendix O The result reads 



n>= 1 + Q2(a) (114) 



where v(a) was defined in Eq. (|97|) . 
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The two-time self-correlation function of the transverse velocity field can also be computed 
from Eqs. (11051) and (11 101) . Again, details of the calculations are given in appendix O For 
s — s' > 0, it is obtained 

V 2 (W + r)Ak 2 

(6u ±ji (k, s)6co ±tj (k', s')) H ~ W*-* x?d eX±mS ~ S ' (U5) 

where 

n OO 

rj x = I dsG{\s\)[e- x ^ k)s -e x ^ k)s ] (116) 
Jo 

It is worth to stress that this result only holds after a transient interval s — s', and for this 
reason it does not reduce to Eq. (I112p for s = s'. On the other hand, the coefficient rj\ can 
be expected to be small, since the second factor in the integrand remains small for the decay 
of the first one. 

To check the theory developed along the paper, Eqs. (11151) and (11121) has been used to 
measure the shear viscosity rj and rf by means of molecular dynamics simulation of dilute 
granular gases. The results turn out to be in qualitative agreement with the theory, in the 
sense that the scaled one-time correlation function is independent of the variable s, and the 
two time correlation function only depends on the difference s — s' and decays exponentially 
after a short transient period. More details of the simulation method employed and the 



analysis of the data is given in 



341 ] . The comparison between the values of rj and rj', 



obtained from the simulation results and the theoretical predictions given by Eqs. (1961) and 
(11141) . respectively is shown in Fig. [TJ Of course, all the simulation data have been obtained 
with low density systems in the HCS. It can be observed that the agreement is very good 
over a quite range of values of the restitution coefficient a, then providing a very strong 
support for both the theory developed here and the specific algorithm used to compute the 
coefficients. 



VII. SUMMARY AND DISCUSSION 

The primary objective here has been to investigate hydrodynamic fluctuations in the 
homogeneous cooling state of dilute granular gases, modeled as an ensemble of inelastic hard 
particles. From this point of view, the fluctuating balance equations (I60l) -(l62l). together with 
the fluctuating Boltzmann equation (14"3"1) provide a solid starting point. The remaining task 
is to construct explicit expressions for the fluctuating flux and the cooling rate in terms of the 
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1.3 




a 



FIG. 1: (Color online) The shear viscosity rj and the new coefficient rf determining the transverse 
velocity fluctuations in granular gases in the HCS. The solid and dashed lines are the theoretical 
predictions for rj and rf given by Eqs. (|96p and (|114p . respectively, normalized by the elastic value 
of the shear viscosity, rj e . The circles (rj/rj e ) and squares (rf /rj e ) are molecular dynamics simulation 
results. The dotted line is the result obtained in the white noise approximation, Eq. (|119p . 

(fluctuating) hydrodynamic fields by using, for instance, the Chapman- Enskog procedure. 
This part of the analysis turns out to be technically rather complex, and has been limited 
here to the particular case of the transverse flow field. 

The structure of the fluctuating balance equations is similar to those for elastic, molecular 
systems with two main differences, both in the equation for the energy, Eq. (IrJiZl) . The 
equation contains a term, 5(o, associated to the fluctuations of the cooling rate and also an 
intrinsic noise term S e . Both give contributions even to zeroth order in the gradients and 
therefore play a relevant role in describing the fluctuations of global properties of the system 

n 

With regards to the fluctuating transverse velocity field, it has been found that it can be 
described by a Langevin equation, but exhibiting two crucial differences as compared with 
the elastic case. The form of the fluctuation dissipation relation changes both qualitatively 
and quantitatively. The second moment of the noise term is not determined by the shear 
viscosity. In addition, the noise is not white, i.e. it presents memory effects. Both aspects 
have been confirmed by the results obtained by molecular dynamics simulations. 
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It could be wondered at what extension the memory effects mentioned above are relevant. 
Suppose the hypothesis of a white noise would have made and Eq. (11051) were substituted 

^(fe, s)W ±J (V, s')) H = 5 id 5 k> _ k ^k 2 7f'S(s - s'), (117) 

with 

POO 

r?' = 2 / dsG(\s\). (118) 
Jo 

By using the same method as outlined in appendix [C] it is found that 



n" = 1 + Q2(a) (119) 

This coefficient is also plotted in Fig. [TJ and it is seen to clearly underestimate the amplitude 
of the second moment of the noise measured in the simulation. It is worth to stress that 
the violation of the elastic fluctuation-dissipation relations is already significant for values 
of the restitution coefficient a of the order of 0.95. 
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APPENDIX A: DERIVATION OF THE FLUCTUATING HYDRODYNAMIC 
EQUATIONS 



In this appendix the derivation of Eqs. ( I60l) - (l62l) will be outlined. The calculations are 
facilitated by using that, for arbitrary functions f(ci, c 2 ) and g(ci, c 2 ), it is 

J dc x J dc2f{c 1 ,C2)T (c 1 ,c 2 )g(c 1 ,c 2 ) 

dc\ \ dc 2 g(c 1 ,C2)T (c 1 ,C2)f(c 1 ,c 2 ), 



where 



T {ci, c 2 ) = I da- 6(ci2 • <r)c 12 ■ a [fe CT (ci, c 2 ) - 1] . 



(Al) 



(A2) 
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The operator bo- is the inverse of b^ 1 defined in Eq. ([6]). From Eqs. (|57|) and (|43j) . 

— 6p(l,s) = / dc—SF^^s) = - dcc- — 6F 1 {l,c,s) + / dcA(c)SF 1 (l,c,s) + S p (l,s), 

(A3) 

with 

= J dcS{l,c,s). (A4) 
By using the definition given in Eq. ( 1281) and the property (1A1I) it is easy to see that 

J rfcA(c)5F 1 (Z,c,s) = 0. (A5) 

Moreover, 

(S p (l, s)S p (V, s'))h = J dci J dc 2 (S(l,c 1 ,s)S(l,c 2 ,s)) H 

= 5{s - s')S{l - V)n- l \- d J d Cl J dc 2 r( Cl , c 2 ), 

(A6) 

where Eqs. (132]) and (153"]) have been employed. But, use again of Eq. ( 1A1I) yields 

da I dc 2 r(a,c 2 ) = o (A7) 



and consequently the second moment of S p vanishes implying that the noise itself also 
identically vanishes. When all the above results are substituted into Eq. (1A3|) and the 
definition in Eq. (|58l) is taken into account, Eq. (|60]) follows directly. 

The balance equation for the fluctuating velocity field is obtained in a similar way. Mul- 
tiplication of Eq. (143]) by a and integration over it gives 



— 6w{h, s) = -— ■ 5U(h, s)+ I da aA(a)SF 1 (h, a, s)S u (h, s), (A8) 

where 11 is defined in Eq. ( 163]) and 

S u {h,s) = J daaS(h,a,s). (A9) 

A simple calculation, using Eq. (lAlh and taking into account the momentum conservation 
in collisions leads to 

J daaA(a)5F l (h,a,s) = ^6u>(l 1 ,s). (Aio) 
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For the correlation of the noise by means of Eq. fl52|) it is found 

(S w (Ji,s)S w (ii,s)>H = -5(s-s')5(h-l' 1 )n- 1 X- d J dc x j dc'^d, {[A( Cl ) + A(t^)] 

x5( Cl - c' 1 )x(c 1 )T (c ll c' 1 )x(ci)x(c' 1 )} . (All) 

The calculation of the first term on the right hand side of the above equation is straightfor- 
ward giving 

J d Cl J dc' lCl c[ [A( Cl ) + A(ci)] (J(ci - ci)x(ci) = Co y dd c lClX (ci) = 1 1, (A12) 

where I is the unit tensor of dimension d. The evaluation of the last term on the right hand 
side of Eq. (IA11I) involves calculating several standard angular integrals. The result is 

J dc x J dc^ctc^Toict, d x )x{ci)x^i) = ^\ T ^d^ d 1 J dc ± j dc! AvX(ci)x(ci) = j I- 

(A13) 

Substitution of Eqs. ( 1A12I) and ( 1A13I) into Eq. ( 1A11I) shows that the second moment of 



the noise vanishes implying that the noise itself also vanish identically. Then Eqs. (1A8I) 
and flATOl) lead to Eq. (|6Tj). 

Finally, multiplication of Eq. fj43|) by 2cf/d and integration over c\ yields 
d . . 2 d .... , d + 2 d . . 

+^ J dc 1 clA(c 1 )SF 1 (l 1 ,c 1 ,s)+S e {h,s), (A14) 

with the heat flux 8(f) defined by Eq. Given that the kinetic energy is not conserved in 
collisions, the calculation of the second term on the right hand side of the above equation is 
more involved than in previous cases, although it is still quite easy to get 

2 f dcic 2 1 A(c 1 )5F 1 (l 1 ,c 1 ,s) = -8Co(l,s) + (o5e(l,s). (A15) 

The expression of the fluctuating cooling rate S( is given in Eq. (1671) . 

Also the calculation of the second moment of the noise term S e is more cumbersome than 
for the other ones. Actually, it is the only one that does not vanish. The result is provided 
by Eq. (J69l) with 



^/ \ 71-2 f r i a /, r \ 2i ^2 i (1 — a 2 )(d + 5 - 2a 



+2(1 + a) (2d + 3 - 3a)c 12 (C • c 12 ) 2 } , (A16) 
where c 12 = c x — c 2 and C = (ci + c 2 )/2. 
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APPENDIX B: DERIVATION OF EQ. flT05|) 



The Fourier transform of Eq. (1521) is 

(S(k, c, s)S(k', c', s '))h = $k,-k>5(s - s')^r(c, c'). (Bl) 

where V = VA~ d , being the volume of the system. Then, for any two components of Wj_, 
from Eq. (I104p it is obtained that 

{W ±>i (k,s)W± jj (k / ,s')) H = ^5 i:j S k ^k 2 J ds x J ds\ j ' dc J dc / c^c ±;i ^c' ±j l((k,c,s 1 ) 

xU(-k, c', 8[)5{8 -s 1 -a' + s' 1 )V±V' A _f(c, c). (B2) 

Here V' is defined like V, but acting on functions of the velocity c', and c\\ denotes de 
component of c along k. Now, it is taken into account that: i) Because calculations are 
restricted to the lowest order in k, which is k 2 , the k factors in the U operators can be 
neglected, i.e. the operator L(k,c) in Eq. (1851) can be replaced by A(c), and ii), Eq. (IHTj) 
implies that, in particular, (x c « c j|£/3) = 0; /3 = 1, . . . , d + 2, for i ^ j. Then, the projection 
operator Vj_ can be eliminated everywhere in Eq. (1B2[) . Note that this can be only done 
after the operator L(k, c) has been replaced by A(c). Then by carrying out the integral over 
Si, it is obtained 

(W ± , l (k,s)W ±J (k / ,s')) H = ^5 k ,-vk 2 J dc J dc ! c\\c ± ^ ± / s - s '^ 

r s ' 

x / dsie"it A W +A ^r(c,C'), (B3) 
Jo 

for s > s'. Again, it is assumed that the HCS is linearly stable with respect to homogeneous 
perturbations for the system under consideration. Then, taking the limit s' ^> 1, the above 
equation can be rewritten as 

(W ± , l (k,s)W ± , 3 (k\s')) u = ^r5 ki „ k ,k 2 J dc J dc'ciic^c'^xe^'^^HCS^c'). (B4) 

The function $hcs( c 5 c is the solution of the equation 

[A(c) + A(c')] ^hcs(c, d) = -f (c, c'), (B5) 

i.e., (see Eq. (E2D), 

#hcs(c,c') = y dlh(l,c;c'). (B6) 
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To give account of the time dependence in Eq. (1B4[) introduce 

iPhcs{c, s; c) = J dlh(l, c, s; c'). 

It obeys the equation (see Eq. (130}) ) 

9 



(B7) 



9s 



A(c) 



^Phcs(c, s; c) = 0. 



(B8) 



By writing Eq. (IB4I) in terms of ipHcs{ c i s \ c ')i Eq. (11051) follows after a trivial change in the 
notation, taking into account the isotropy of A(c). 



APPENDIX C: SOME DETAILS OF THE CALCULATIONS DESCRIBED IN 

sec. ED 

The explicit expression of the function G(|s|) defined in Eq. (11061) is 

G(\s\) = J dc J dc'c x c y c' x c' y e sA{c) x(c)S(c-c') = J dcc x c y e sA{c) X (c)c x c y . (CI) 
The term involving the two-particle correlation function gucs in Eq. (11091) has been ne- 



33] has 



glected. The reason for this is that numerical simulations of dilute granular gases 
shown that it gives much smaller contributions than the term kept. Substitution of the 
above expression in Eq. ( 11131) and integration over time s leads to 



i] — 2 J dcc x c y C X y(c) , 



(C2) 



where C xy (c) obeys the equation 



A(c) + I' 



C xy {c) = -x(c)c x c y . 



(C3) 



This equation will be now solved in the first Sonine approximation, i.e. by keeping only the 
first term of the expansion of the solution in Sonine polynomials 
of the operator A(c) has the form 



20[, that due the isotropy 



C X y(c) = AC X Cylf (c) 



(C4) 



where (pa is the Gaussian 



Vo = 1T - d / 2 e- c2 . 



(C5) 
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To determine the constant A, Eq. (1C3j) is multiplied by c x c y and integrated over c, after 
replacing C xy (c) by its expression in Eq. (IC4I) . Moreover, the velocity distribution of the 
HCS is also approximated by its expression in the first Sonine approximation, Eq. (1221) . 
From this point, the only remaining task is to evaluate several integrals. Although lengthy, 
the calculation is straightforward and similar to the one needed to compute the Navier-Stokes 
transport coefficients 30j. The result found for the constant A reads 

1 + a 2 (a) 



A 



(C6) 



with v{a) given by Eq. ( 1971) . Finally, Eq. ( 1C4I) is substituted in Eq. (1C2I) and the result 
( 11141) follows quite easily. 

In order to derive an expression for the two-time correlation function of the transverse 
velocity, Eqs. (11051) and (jllOp are employed. Taking s > s' 1, 



V 2 

(6u± !i (k, s)5u± d (k' ', s')) H = -j^- k 2 5 itj 5 k - k > 



d Sl / ds[ e x ^ 2s '- Sl ~ s ^G(\ Sl - 41) 



+ 



d s ' ie Ax(fc)(»+«'-»i-«i) G (| Sl _ s ' i |) 



(C7) 



The first integral on the right hand side is related with the coefficient v( computed before, 



oo J — oo 



d Sl / ds; e^W^-^-^Gflsi - si|) 



(C8) 



2A ± (ifc)' 

The other integral in Eq. ( 1C7I) can be evaluated for instance by introducing the variables 
g = si — Sj and S 1 = (si + s' l )/2. Then, the integral is decomposed into three. For large 
enough s — s', the only remaining one is 



dg G(\g\) 



dS e X±(k)(s+s'-2S) 



Jl e \±{k)(s-s>) 



io Js'-i 2X±(k) 

with ft given by Eq. (fTT6l) . Use of Eqs. (JC8]) and (JC9} into Eq. (jDZj) gives Eq. ( flT5j) . 



(C9) 
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